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Abstract 


We report results (on an intermediate statistics sample) of a study of weak 
semi-leptonic formfactors of B and D decays, addressing the uncertainties 
from mass extrapolations to chiral and to heavy quarks. Moreover, we present 
a nonperturbative test to the LMK current renormalization scheme for vector 
current transition matrix elements and hnd remarkable agreement. 


I. INTRODUCTION 

The determination of weak formfactors for the semi-leptonic decays of heavy-light mesons 
presents a big challenge for further exploration of the flavordynamics within the Standard 
Model, both experimentally and theoretically. The reason is that the study of weak in¬ 
teraction amplitudes like the Kobayashi-Maskawa parameters calls for good control of the 
hadronic binding effects in the analysis. 

Heavy quark effective theory provides a suitable framework of understanding the ’ge¬ 
ometric’ part of the spin-flavor structure of heavy quark decay in the infinite mass limit 


Hnvited talk given at the German-Japanese workshop Physics on Parallel Computers, Yamagata 
1995 

*Work supported by DFG grant Schi 257/3-2 and 3-3 


1 



Mh —^ oo, but cannot predict the l/M^ corrections in detail which are likely to be impor¬ 
tant in the mass regions of D- and B-mesons. For this reason great effort has been spent 
ever since the pioneering work of the El-Khadra et al [1] to compute QCD effects on these 
decays by means of lattice gauge theory. 

In this contribution, we present a status report about our ongoing work on 24^ X 64 
lattices with quenched Wilson fermions at /3 = 6.3, based on 60 conhgurations (out of a 
total sample of 100 conhgurations) that we are carrying out on the 32 node connection 
machine CM5 with 40 Gbytes of parallel disc space SDA, at Wuppertal University. The 
light quarks were represented by the following set of hopping parameters nf. .1511, .1507, 
.1490, .1450, while the heavy quarks were varied over the set of Hh\ .1400, .1350, .1300, and 
.1200. The initial meson was chosen at rest, while the hnal meson could carry 11 different 
momenta (in lattice units): (0,0,0), (1,0,0), (1,1,0), (1,1,1), (2,0,0), and permutations. 


II. MASS EXTRAPOLATIONS 

Let us for the moment consider the momentum transfer dependence of the formfactors, 
leaving aside the issue of the lattice renormalization of these quantities. 

Chiral extrapolation Fig.l gives an idea on the quality of our data, in terms of the 
conventional formfactors /_|_ and /o for the transition PS —> PS + leptons as well as the 
formfactors Ai, A 2 , and V for the decay PS —)• U + leptons. For the formfactor dehnitions, 
see Ref. [3]. Kinematically, this hgure refers to a heavy quark (kj- = .1350 decaying into a 
light one [k, = .1490), with a spectator quark, = .1490. The full lines represent the 
pole dominance predictions, with pole masses from the lattice two-point functions. 

In Fig. 2 we demonstrate, for comparison the effects of chiral extrapolation for kj and 
Kspec- Physically this hgure relates to the decays D ^ p and U —)• tt. 

Unfortunately, our lattice allows only for a small variation in the momentum transfer 
between the initial and hnal state hadrons, q^. Generally, we hnd consistency with (but no 
compelling evidence for) pole dominance. 

Extrapolation to B Meson It is well known that currently accessible lattice sizes 
do not allow for a direct simulation of B decays and therefore require a substantial mass 
extrapolation, from the D to the B meson regime. While the transitions i? —)• U, D* lend 
themselves rather easily to a heavy quark effective theory (HQET) like analysis [7], the 
situation is much more involved in the decays B ^ tt, p. One option is to carry out this 
extrapolation in a manner suggested by HQET, i.e. at hxed value of cu = ViVj [3,8]^. This 
extrapolation implies a big change in which has to be compensated by a substantial 
counter-extrapolation to q^ = 0, the conventional point of comparison. Given the uncertain¬ 
ties in the formfactors this procedure introduces sizeable systematic errors. We therefore 
prefer to carry out the extrapolations at fixed = 0, using different functional forms of 
Mps dependence. This appears to be a better controlled way to judge the systematic errors 
on the hnal numbers. 

In Fig.3 we display four of the formfactors F{q^ = 0) versus the inverse mass of the 
decaying meson, Mfig. The hnal meson corresponds to K and JF*, respectively. The differ- 


^In ref. [8] also an extrapolation in mass at = 0 has been performed. 
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ent curves represent different types of linear extrapolations to the B meson mass (leftmost 
points): full lines fits to F itself, dashed lines hts to F\jMps^ dotted lines hts to F\jMps ^■ 
The plots illustrate that the data is not sensitive yet to discriminate between the M depen¬ 
dencies of the ansatze. 


III. RENORMALIZATION 

In order to relate our results to the continuum we need to renormalize the matrix elements 
of the local (non-conserved) vector and axial vector currents. It has been repeatedly pointed 
out that - on top of their short distance part - the renormalization factors may receive large 
0{ma) contributions at large ma. 

An obvious starting point is to take these contributions into account by approximately 
adjusting the normalization of the lattice quark propagator A(l,p = 0) to the one of the 
continuum quark propagator [2] which has been advocated and rehned more recently by 
Lepage, Kronfeld and Mackenzie [5], and will hence be quoted as LMK scheme. From the 
free case relation A(l,p = 0)“*^* = 2k{1 + ma)A{t,p = 0)^“** one uses mean held arguments 
to obtain 

A{t,p = 0^ ~ 2KNLMK{K)A{t,p = 0)'“" (1) 

for the interacting case. As pointed out in ref. [6]. the normalization factor = 

^(1 — 1^) can be “eaten up” by an additional time hop of the lattice quark propagator 

A(l,p = ~ A(l +l,p = 0)'““ . (2) 

In this section we will measure the mass dependence of the renormalization constant of 
the time component of the vector current and compare the results to the LMK predictions 
as well as to the naive expectations, in which the LMK factor is absent. 

In order to improve on statistical accuracy, we will work with ratios of 3-point functions, 
rather than those functions themselves. We emphasize that the mass dependence of such 
ratios directly measures lattice artefacts, as it disappears in the continuum limit. 

We consider the ratio^ 

„ _ {PS{mHp mi))p 

^ {PS{mHprni)\V^°‘^\PS{mH,rni))p 

The ratio Rhh{^) yields directly the renormalization factor of as is unaltered 

in this instance. For the study of decays we are rather interested in the ratios i?/i//(0) and 
Rhh{p)^ however. 

In order to derive predictions for the mass dependence of the ratios R above we have to 
study the 3-point functions 

2/0, Xq) = 0) ^ yo)PHl{x, Xq)) , (4) 


^A similar analysis has been carried out in ref. [9]. 


(3) 
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where P denotes the operator of the pseudoscalar meson and 


yo) = yoho'^’niy, yo) (5) 

yo) = ^ l^H'iy, yo)ipo - l)Uoiy, ?/o)^ff(y, yo + 1) 

+ ^H'{ypjo + l){7o + l)U^{ypjo)'i’H{ypjo)] ■ ( 6 ) 

For 0 ^ j/o ^ 3;o one hnds Rhh' = G^hh'iK^^hh'I- terms of quark propagators, A, one 
easily obtains 

GHH,i{0,yo,xo) = (^Tr {A/j(y,j/o,®, a:o)75A/(F,xo, 0,0)A];^,(y,j/o, 0 , 0 ) 7570 }) (7) 


?/0, Xq) — (- X] [ (®) 

x,y 

Tr {A//(y, yo + 1, F, 2 : 0 ) 75 A;(F, xq, 0, 0)A];^,(y, j/o, 0, 0 ) 75(70 - l)Uo{y, J/o)} 

+ Tr {A//(y,j/o,®, 2:0)75 A/(F, 2:0,0, 0 )A];^,(y,j/o + 1,0,0)75(70 + l)U^{y,yo)}]) ■ 


Inserting the relations (1) and (2) for the quark propagators one can now read off the 
LMK predictions for the mass dependence of the ratios i?(0). Note that the second term in 
eq.( 8 ) is expected to be negligible within our mass range^. Because of j/o ~ 2:0 < 0 one gets 
AHiyo,xo) (1 — 70 ) for very heavy quarks, and the second term becomes proportional to 
(1 — 7o)(l + 7 o) = 0. The LMK predictions for the ratios therefore read 


Rhh{^) 




V 


3 

4 Kc 


2Kff 


RhH(O) 


yloc 


ycons 


3 

4 Kc 


2kh 


(9) 


where Zy‘^ and Zy”® are constants. We emphasize that Rhni^) only depends on niH but 
not on rrih- This is a consequence of the cancellation of the second term in eq.(8). The naive 
prediction simply reads Rhh{^) = = RhH- 

In hg.4 we compare our numerical results for Rhh{^) with the LMK and the naive 

predictions. We have inserted the (1-loop) perturbative values Zy = 1 — 0.82/47r5'^, Zy‘^ = 
1 — 0.1740815'^, using the boosted coupling and assumed Zy”® = Zy”® = 1. Throughout 
our mass range of niH we hud the LMK prediction to approximate our data within 5%. The 
naive prediction, however, fails drastically. 

For actual physics applications it is more interesting to study the situation, where initial 
and hnal nonspectator quarks carry different masses. For that matter, hg.4 shows the ratio 
7?/i_ff(0) as a function of nih at hxed niH. Again we hud good agreement with the LMK 
prediction. This hnding has important bearing on the reliability of our calculations with 
light hnal states. 


®This has also been verified numerically from our data. 
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There remains the hnal issue, whether LMK holds as well in the case of unequal momenta 
of initial and hnal state hadrons. To this end, hg.4 displays the ratio Rhh{p) on the 
momentum set \pa\ = 27r/64 X {0,l,\/2,2}- The solid line represents the above LMK 
prediction, for p = 0. We hnd that the data moves only marginally with increasing p and 
is therefore well approximated by the LMK prediction over the entire p range. This result 
is not at all trivial, as one would expect the normalization factor to be affected by the 
momentum. 


IV. PHYSICS RESULTS 

In table I we have collected our results for the B decays into charmed particles, while 
the tables III and II contain our results on decays into light hadrons and on D decays, 
respectively, together with the data of other lattice investigations, as well as model and sum 
rule calculations. 

We hnd good agreement with experimental data for the D —)• JL, K* decays, with possible 
exception of ^ 2 ( 0 ), where the errors are large. Note that the kinematical constraint /+(0) = 
/o(0) is well fulhlled by our data! 

While our results for D ^ tt, p are generally fairly consistent with the results from 
other lattice investigations, there is still rather large scatter in the results for the decays 
B ^ TT, p. This is related to the heavy mass extrapolations. We have not imposed the 
constraint /+(0) = /o(0) in our analysis. Therefore the difference in these quantities rehects 
systematic uncertainties. The extrapolative situation for these formfactors can be traced 
from the upper plots in hg.3. /o is determined with higher precision, due to the kinematical 
point pf = 0. A ht using the constraint /+(0) = /o(0) along the mass trajectory would 
favour the value 0.20, of course. A similar kind of uncertainty is borne out by the range 
quoted recently by the UKQCD collaboration [8] for /+(0) = /o(0). They have used the 
latter constraint within an HQET inspired extrapolation procedure, using several variants 
for the (j'^-dependence of the formfactors. 

At this stage, one should note, that reliable knowledge obout the Z-factors is urgently 
needed, as they enter directly into the lattice predictions for the formfactors. It is therefore 
of great interest to continue studies about nonperturbative lattice renormalization. 
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FIGURES 



FIG. 1. dependence of the formfactors, at Ki = .1350, and Kf = K^pec = .1490. The curves 
refer to hts: pole dominance with mass from two-point function (full line), with mass htted to 
formfactor (dashed line), and linear ht (dotted line). 
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1/Mp3 1/Mp3 

FIG. 3. Extrapolation of formfactors in the mass of the decaying pseudoscalar particle [in 

GeV]^ with hnal state hadron K and ItG, respectively. The data are included as crosses, and the 
extrapolated points are plotted at the inverse B-mass, Mpg = .2. All hts are linear: (a) to F itself 
(full line), (b) to FMp^g (dashed line), (c) to FMpg^^ (dotted lines). The interpolated results at 
the D-mass are included for convenience (open symbols). 
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FIG. 4. The ratios Rhh{^)^ RhH{^) and Rhh{p) as a function of the quarkmass. Solid lines 


represent the LMK and dashed lines the naive prediction. 
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TABLES 


TABLE 1. Formfactors for B D, D* und D*. 


B -r D,D* 


/+(0) 

/o(0) 

E(0) 

Ai(0) 

A2(0) 

1.01(60) 

0.88(29) 

1.84(1.10) 

0.97(32) 

2.76(1.58) 


B,^D„ D* 


/+(0) 

/o(0) 

E(0) 

Ai(0) 

A2(0) 

0.99(42) 

0.85(24) 

1.84(1.10) 

0.91(24) 

2.25(1.12) 
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TABLE II. Formfactors at = 0 for the decays D K und D K*. 



Method 

/+(0) 

/o(0) 

E(0) 

Ai(0) 

A2(0) 

World av. [10] 

Exp’t 

0.77(4) 


1.16(16) 

0.61(5) 

0.45(9) 

This work 


0.71(12)t].° 

0.67(6)t? 

1.34(24)tJ® 

0.61(6)t? 

0.83(20)^7 

LMMS [11] 

Lsjtt . 

0.63(8) 


0 .86(10) 

0.53(3) 

0.19(21) 

BKS [12] 

Wilson 

0.90t®t2i 

0 7n+^+24 

U. / U_g_24 

1 /LQ+45-1-48 

J-.4:0_45_4g 

n oq + 14 + 28 
U.00_;L4-28 

n 50 + 14-1-24 
U.O»_44_23 

ELC [3] 


0.60(15)(7) 


0.86(24) 

0.64(16) 

0.40(28)(4) 

BG [13] 


0.63(3) 

0 .66(2) 

1.18(7) 

0.65(3) 

0.41(7) 

APE [14] 

Lsjtt . 

0.72(9) 


1 .0(2) 

0.64(11) 

0.46(34) 

UKQCD [15] 

Clover 

0.67^1 

0.65(7) 

1 .0lt?° 

0.70t7 

0.66^7 

ISGW [16] 


0.76 - 0.82 


1.1 

0.8 

0.8 

WSB [17] 

Quark- 

0.76 


1.27 

0.88 

1.15 

KS [18] 

Models 

0.76 


0.8 

0.82 

0.8 

GS [19] 


0.69 


1.5 

0.73 

0.55 

BBD [20] 


0.60tJ® 


1.10(25) 

0.50(15) 

0.60(15) 

AEK [21] 

Sumrules 

0.60(15) 





DP [21] 


0.75(5) 
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TABLE III. Formfactors at = 0 for the decays D and i? to vr and p. 




D 

^ TT und D - 

P 




Method 

/+(0) 

/o(0) 

U(0) 

Ai(0) 

7l2(0) 

This work 


0.68(13)1].° 

0.65(5)1? 

1.31(25)1]? 

0.59(7)1? 

0.83(20)1?^ 

LMMS [11] 

Lsjtt. 

0.58(9) 


0.77(9) 

0.47(7) 

-0.07(42) 

BKS [12] 

Wilson 

0.84(12)(35) 

0.70(8)(24) 

1.07(49)(35) 

0.65(15)1?? 

0.591?? 

BG [13] 


0.52(6) 

0.57(4) 

1.12(11) 

0.58(4) 

0.34(10) 

UKQCD [15] 

Clover 

0.6ll]J 


0.951]? 

0.631? 

0.511?? 

ISGW [16] 


0.51 

0.51 




WSB [17] 

Quark- 

0.69 


1.23 

0.78 

0.92 

KS [18] 

Models 

0.69 


1.23 

0.78 

0.92 

AEK, DP [21] 

Sumrules 

0.6-0.75 

0.6-0.75 






B 

^ TT und B - 

^ P 




Method 

/+(0) 

/o(0) 

U(0) 

Ai(0) 

4l2(0) 

This work 

Lsjtt. 

0.50(14)1^ 

0.20(3)1? 

0.61(23)1? 

0.16(4)1?? 

0.72(35)1?° 

ELC a [3] 

Wilson 

0.28(14) 


0.37(14) 

0.24(6) 

0.39(24) 

ELC b [3] 


0.33(17) 


0.40(16) 

0.21(5) 

0.47(28) 

APE a [14] 

Lsjtt. 

0.29(6) 


0.45(22) 

0.29(16) 

0.24(56) 

APE b [14] 

Clover 

0.35(8) 


0.53(31) 

0.24(12) 

0.27(80) 

UKQCD [8] 


0.24tt 

0.241? 


0.271? ± 3 


ISGW [16] 

Quark- 

0.09 


0.27 

0.05 

0.02 

WSB [17] 

Models 

0.33 


0.33 

0.28 

0.28 

Ball [22] 

Sumrules 

0.26(2) 


0.6(2) 

0.5(1) 

0.4(2) 
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